Minimal set of local measurements and classical communication for two-mode 
Gaussian state entanglement quantification 
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We develop the minimal requirements for the complete entanglement quantification of an ar- 
bitrary two-mode bipartite Gaussian state via local measurements and a classical communication 
channel. The minimal set of measurements is presented as a reconstruction protocol of local covari- 
ance matrices and no previous knowledge of the state is required but its Gaussian character. The 
protocol becomes very simple mostly when dealing with Gaussian states transformed to its stan- 
dard form, since photocounting/intensity measurements define the whole set of entangled states. In 
addition, conditioned on some prior information, the protocol is also useful for a complete global 
state reconstruction. 
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Quantum communication protocols extend the infor- 
mation theoretical notion of channel [l[ to the quan- 
tum domain by incorporating non-local entangled states. 
Those channels are generated by the preparation of a 
pair (or more) of quantum systems in an entangled state, 
which are then separated to establish non-local correla- 
tions , allowing several communication tasks otherwise 
unattainable via classical channels However, for most 
of the quantum protocols to work properly (deterministi- 
cally) one has first to be able to prepare maximally pure 
entangled states and then to guarantee that those states 
stay pure or nearly pure during all the processing time. 
An important problem then arises in this whole process: 
One has to check the "quality" (the amount of entangle- 
ment and purity) of the quantum channel, while usually 
the only available tools for that are local measurements 
(operations) and one (or several) classical channel. 

The quest for an optimal and general solution for this 
problem has generated a vast literature on the character- 
ization of entangled states under local operations and 
classical communication (LOCC), either for qubits Q 
or for continuous variable systems of the Gaussian type 
[^, @|. Gaussian states (completely described by up to 
second order moments) are particularly important since 
they can be easily generated with radiation field modes. 
Moreover, operations that keep the Gaussian character 
(so-called Gaussian operations) are given by the transfor- 
mations induced by linear (active and passive) optical de- 
vices (beam-splitters, phase-shifters, and squeezers) Q. 
A particular result for this kind of state is that it is im- 
possible to distill entanglement out of a set of Gaussian 
states through Gaussian operations 

Assuming one is left with only Gaussian local opera- 
tions and a classical channel (GLOCC), how is it possible 
to infer the quality of a quantum channel in use? For a 
two-mode Gaussian state one possibility is to access di- 
rectly the entanglement properties of the system after a 
proper manipulation of the two modes [^, 0] . This proce- 
dure requires, however, that the two parties (modes) be 
recombined in a beam-splitter (non-local unitary opera- 



tion) in which their entanglement content are transferred 
to local properties of one of the output modes. Another 
possible way is to completely reconstruct the bipartite 
quantum system, a resource demanding task [lo| which 
also requires global operations here forbidden. 

In this Letter we demonstrate a minimal set of GLOCC 
to completely quantify the entanglement of a two-mode 
Gaussian state. As a bonus of this procedure one can 
also assess the purity of the Gaussian state and, for some 
particular classes of states, reconstruct the bipartite co- 
variance matrix. The protocol consists mainly in the at- 
tainment, via local measurements, of all the symplectic 
invariants that allows, for example, one to test the sep- 
arability of the system, to know its P-representability 
properties, and to quantify its entanglement content. We 
also show that for a particular class of Gaussian states 
belonging to the set of symmetric Gaussian states [11 



the Einstein-Podolsky-Roscn (EPR) states and general 
mixed squeezed states , the protocol becomes straightfor- 
ward due to the relative easiness one obtains the correla- 
tion matrix elements from local measurement outcomes. 
Moreover, since P-representability and separability for 
these kind of states are equivalent, we show that for two- 
mode thermal squeezed states with internal noise 12] it 
is possible to decide whether or not they are separable 
via local photon number measurements. 

A two-mode Gaussian state 
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is characterized by 
its Gaussian characteristic function C{cx.) = e"^"*'^"^ 
where = (a*, ai, 02) are complex numbers and 
fli (a|) and 02 (aj) the annihilation (creation) opera- 
tors for parties 1 and 2, respectively [13[. The covari- 
ance matrix V describing all the second order moments 
V^j = {-iy+^ {v^v] + v]vij/2, where (i;i, 1^2, '^3, 1^4) = 
(ai, a|, 02, ftj), is given by 



V 



Vi c 
Ct V2 



/ 


ni + \ 


mi 


rris 










ni + i 


m* 


m* 






m* 


rric 


"2 + i 


m.2 




V 


m* 


rUs 


m* 


n2 + 





2 



Vi and V2 are local Hcrmitian matrices while C is the 
correlation between the two parties. Any covariance ma- 
trix must be positive semidefinite V > and the gen- 
eralized uncertainty principle, V -t- (1/2)E > 0, where 
E = diag(Z, Z) and Z = diag(l, -1), must hold 

From local measurements on both modes of pi2, ei- 
ther through homodyne detection (see [l^l and references 
therein) or alternatively by employing single-photon de- 
tectors [1^, the local covariance matrices Vi and V2 
can be reconstructed. Remark that for the reconstruc- 
tion of the global matrix V, and therefore the joint bi- 
partite state, one has to obtain C. Obviously, global joint 
measurements achieved through recombination of the two 
parties in a beam-splitter followed by local homodyne de- 
tections are forbidden. Thus one has to deal only with lo- 
cal measurements whose results can be sent through clas- 
sical communication channels to the other party. As we 
now show, there are minimal operations/measurements 
that can be performed locally on the system to attain 
I det C| and detV. These quantities, together with det Vi 
and det V2 , will be shown to be all that one needs to de- 
termine whether or not a two-mode Gaussian state is 
entangled as well as how much it is entangled. As it will 
become clear, the required set of operations is minimal 
in the sense that only two local measurement procedures 
are needed - one to characterize local covariance matri- 
ces and another to locally assess the parity of one of the 
modes. 

First of all let us introduce an important result [l^ . 
Given a two-mode Gaussian state with density operator 
P12 and covariance matrix V we can define the Gaus- 
sian operator cti = Tr2 |e''^°2'^2pi2| , whose covariance 
matrix Fi is the Schur complement [3| of V relative to 
V2: 

Ti^Vi-CV^^CK (1) 

The meaning of ai is best appreciated through a par- 
tial trace in the Fock basis: cti = Snei.e,i2("l/^i2l'^)2 — 
Sn„di2('^|pi2|n)2 = Pie - Pio: being equal to the dif- 
ference between Alice's mode states conditioned, respcc- 
tively, to even and odd parity measurement results by 
Bob [l7[. While pi^ and pi^ are not generally Gaussian, 
CTi is a Gaussian operator, and Fi can be built with only 
second order moments of these conditioned states. 

Now suppose that Alice and Bob share many copies 
of a two-mode Gaussian state. The protocol works as 
follows: (z) Firstly, in a subensemble of the copies, each 
party performs a set of local measurements in such a 
manner to obtain the covariance matrices Vi and V2, 
corresponding to the reduced operators pi = Tr2{pi2} 
and p2 = Tri{pi2}; (ii) Then Bob informs Alice, via 
a classical communication channel, the matrix elements 
of V2; {Hi) After that, for the remaining copies. Bob 
performs parity measurements on his mode, letting Alice 
know to which copies does that operation correspond and 



the respective outcomes, i.e. even parity (eigenvalue 1) or 
odd parity (eigenvalue -1); (iv) Alice then separates her 
copies in two groups, the even (e) and the odd (o) ones. 
The first group (e) contains all the copies conditioned on 
an even parity measurement on Bob's copies. The other 
one (o) contains all the remaining copies, namely those 
conditioned on an odd parity measurement at Bob's; {v) 
For each group, Alice measures the respective correlation 
matrices Vie and Vio; (vi) Finally, she obtains Fi (Eq. 
([1])) subtracting the odd correlation matrix from the even 
one [l3|: Ti = Vie - Viq. Remarkably, with Vi, V2 
and Fi in hand Alice is able to completely characterize 
the Gaussian state's entanglement content as well as its 
purity without any global or non-local measurements. 

Remembering that a two-mode Gaussian state's purity 
V is equal to l/(4\/det V) ^ and using the identity ^ 

detV = dctV2detFi, (2) 

Alice readily obtains the purity of the channel: V = 
l/(4Vdet V2 det Fi). Her next task is to decide whether 
or not she deals with an entangled two-mode Gaussian 
state. Using the Simon separability test she knows 
that it is not entangled if, and only if, 

/1/2 + (1/4 - I/3I)' -h> ih + /2)/4, (3) 

where Ii = det Vi , I2 = det V2 , I3 = det C , and 
/4 = tr(ViZCZV2ZCtZ). These four quantities are the 
local symplectic invariants, belonging to the Sp{2, R) (E) 
Sp(2,R) group [20I, that characterizes all the entangle- 
ment properties of a two-mode Gaussian state. Alice 
already has Ii and l2- We must show, however, how she 
can obtain I/3I and I4. Since one can prove that Q 

h^2\l3\^/hh, (4) 

we just need to show how I/3I is obtained from /i, I2, 
and Iv — dot V, the three pieces of information locally 
available to Alice. To achieve this goal we first note that 
a direct calculation gives ly = I1I2 ^ I4 + li- Using 
Eq. dl]) we see that I/3I follows from I/3P - 2I/3IV/1/2 + 
/1/2 — /y = 0. One of its roots is not acceptable since it 
implies V < 0. Therefore, we are left with 

\h\ = Vhh- Vh^- (5) 

Hence, substituting Eqs. (|4|) and ^ in Eq. Alice is 
able to unequivocally tell whether or not she shares an 
entangled two-mode Gaussian state with Bob. Finally, if 
her state is entangled then I3 < [2^ and, for a symmet- 
ric state (/i = 72), Alice can quantify its entanglement 
via the entanglement of formation (Ef) [2ll.[2^: 

Efipu) = f (^2^Ii + \h\~Vh + 2h\h\^ , (6) 

where f{x) = c+(a;) log2(c+(a;)) — c_ (x) log2(c_ (a;)) and 
c±{x) = {x~^l'^ ± x^l'^Y' l\. For arbitrary two-mode 
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Gaussian states (L ^ I2) Alice can work with lower 
bounds for Ej [2l| or calculate its negativity or log- 
arithmic negativity [2^. This last two quantities are 
the best entanglement quantifiers for non-symmetric two- 
mode Gaussian states and are given as analytical func- 
tions [3, 24 1 of the four invariants here obtained from 
local measurements: Ji, /2, I/3I 
h 



h, h, \h\ = Vhh - VTv, and 
2|/3|/rr/2, with Iv = detV given by Eq. It is 



worth mentioning that Ii {I2) can easily be determined 
by the measurement of the purity (Wigner function at 
the origin of the phase space) of Alice's (Bob's) mode 



alone la, l25| . This measurement is less demanding than 
the ones required to reconstruct Vi and V2 0]. 

Besides furnishing all the entanglement properties of 
an arbitrary two-mode Gaussian state, the previous local 
protocol can also be employed to reconstruct the covari- 
ance matrix for some particular types of Gaussian states. 
To see this, let Fi be explicitly written as 
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where 



rji = {a\ai)e - {a\ai)o, (8) 
Ml = (a?)e - (a?)„, Hi = ((al)2), - ((al)^)^, (9) 

being and (•)o the mean values for Alice's even and 
odd subensembles, respectively. From this identity it is 
clear that Fi does not necessarily represent a physical 
state since -qi can take negative values [13] ■ From Eq. 
(fT|) we obtain the following two relations, 
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Eqs. (|10[) and (|11[) give the matrix elements of Fi as a 
function of the matrix elements of V. If nic and ms arc 
real (if either rric or nis is zero) Eqs. pUj) and pT|) can 
be inverted to give uic and ms (cither ms or mc). 

Let us explicitly solve the previous equations for an 
important case, namely the ones in which CC^ = |mipl, 
where i = c or s and I is the identity matrix. The 
states comprehending this class are the ones where C 
has only diagonal or non diagonal elements, i.e., ms ~ 
and TTic ^ or TOc = and ^ 0, reducing the un- 
known quantities to two, namely the absolute value and 
the phase of ms or mc- Remark that if i = s the system 
is separable, since det C = |tosP > 0, i.e., the correlation 
between the two modes is strictly classical Other- 
wise, if i = c the state is not necessarily separable, possi- 
bly being entangled, for in this case det C = — |mcp < 0. 



This last case is more interesting since it represents a 
class of states that might show non-local features [2^. 
From Eqs. and pT|) the diagonal (off-diagonal) 



elements of C, 



|TOi|e"^', for i ^ s {i = c), 
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where m2c = nij and m2s = 7712. Note that whenever 
7ii2 = 0, (jji becomes undetermined. This problem can be 
solved by locally (unitary) transforming the two-mode 
squeezed state to a matrix V2 with 7712 ^ 0, where 4)[ 
can be determined. Then, transforming back, we get 0^. 
Fortunately, there are various experimentally available 
bipartite Gaussian states in which all the parameters are 
real, ms{mc) = mi = 1112 = 0, and mc{ms) ^ 0. For 
these states, Eq. ([T2|) is sufficient to determine C. 

A natural and important example belonging to this 
class is the two-mode thermal squeezed state [12|, which 
is generated in a nonlinear crystal with internal noise. Its 
covariance matrix is 
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where 71 and mc arc time dependent functions having as 
parameters the relaxation constant of the bath as well 
as the nonlinearity of the crystal [l2|. In this case the 
protocol involves only simple local measurements, i.e., 
those to get 71, (a|ai)e and (a\ai)o (or equivalently 771) 
by Alice, and the parity measurements by Bob. The 
classical communication corresponds to Bob informing 
Alice the instances he performs the parity measurement 
in his mode and the respective outcomes. Hence, Eq. p2|) 
reduces to 

r7i2 = (n-7?i)(7i+l/2). (15) 

Experimentally, 71 and 7/1 (Eq. ([8])) are readily obtained 
by photodetection, while the parity measurement is re- 
lated to the determination of Bob's mode Wigner func- 
tion at the origin of the phase-space [2^ , or alternatively 
to his mode's purity, both of which can be measured by 
photocounting experiments 



We can also study the P-representability [27| for the 
state (dH), which in this case is equivalent to the Simon 
separability test [ll|, [2^ . A two- mode Gaussian state is 
P-rcprcsentable iff V— > 0, where I is the unity matrix 
of dimension 4. Explicitly, this separability condition in 
terms of the elements of ([T4|) is equivalent to n > \mc\. 
From this inequality and Eq. (|15p we see that for a given 
71 there exists a bound for 771 below which the states are 
entangled (upper solid curve in Fig. [1]): 

71/2 n/2 



1/2 



< i;i < 



1/2- 



(16) 
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The left bound in Eq. ([16]) (lower solid curve in Fig. 1) 
is a consequence of the uncertainty principle, delimiting 
the set of all physical symmetric Gaussian states (SGS). 
This bound is marked by all the pure states and the up- 
per curve bounds (from below) the subset of all sepa- 
rable (P-representable) states Thus, for the SGS 
class, photon number measurements, before and after 
Bob's parity measurements, are all Alice needs to dis- 
cover whether or not her mode is entangled with Bob's. 
The exquisite symmetry of those two antagonistic bounds 
is quite surprising, and possibly valid only for the SGS 
class. There is another interesting feature for the SGS set 
that should be emphasized. Note that ryi = contains all 
the states where Bob has equal chances of getting even 
or odd outcomes for his parity measurements, delimiting 
two subsets (even and odd). The even subset contains 
all the states where Bob has greater probabilities of get- 
ting even outcomes while the odd subset contains all the 
states where he has greater probabilities of getting odd 
outcomes. The entanglement for states belonging to the 
SGS can be quantified through Ef (Eq. ([6])) as depicted 
by the color scale in Fig. [T] It is remarkable that the 
most entangled states (including the pure ones) are con- 
centrated in the 771 < odd subset. 




FIG. 1: (Color online) Above the upper solid curve lie the 
separable states. Below it, entanglement is quantified via Ef 
(Eq. (|6])) up to the lower curve, where the pure entangled 
states are located. Below this curve there exist no physical 
states. 

In conclusion, we have presented the minimal set of 
local operations and classical communication that allows 
one to quantify the entanglement of an arbitrary two- 
mode Gaussian state. One important step towards the 
derivation of this protocol was the mathematical identity 
relating the two-mode covariance matrix determinant to 
the product of two local quantities, namely the determi- 
nants of the one-mode correlation matrix and its Schur 
complement. In addition, we have also shown that the 
Schur complement of one of the modes' covariance ma- 
trix is obtained via a set of parity measurements on the 
other one. We have also explicitly discussed how the 



protocol works for a particular class of Gaussian states 
belonging to the SGS set. Within this class, for states 
written in its standard form, we have shown that only 
photon number measurements (made before and after a 
parity measurement on the other mode) are needed to 
completely characterize the state's entanglement. 
This work is supported by FAPESP and CNPq. 
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ERRATUM 

Eq. (4), on page 2, of [Phys. Rev. Lett. 98, 150501 
(2007)] is not so general as we have previously thought. 
However, our scheme does not rely on it, as we show in 
what follows. 

The experimental proposal we presented in our Letter 
allows one to locally obtain the matrices Vi,V2, and 
Fi, without assuming any particular form for the co- 
variance matrix V. We now show that with these three 
matrices wc can determine the four invariants that com- 
pletely characterize the entanglement content of a two- 
mode Gaussian state. The first two invariants are 

/i=det(Vi), /2=det(V2). (17) 

The third one is calculated remembering that Fi = Vi — 
CV2 """Ct. A simple algebra on the previous expression 
gives, dct(Vi-ri) = det(C)det(V2^)det(Cl'). But 
det(C) = det(Ct) = /g and det(V2^) = l/det(V2) = 
I//2. Hence 

1/3 1 - ^hdet (Vi-Ti). (18) 
Furthermore, Fi satisfies another mathematical inden- 



tity, Jy = detV = detV2detri. Therefore, since we have 
Fi and V2, wc can also obtain Jy. But /y is related 
to the other four invariants by the following expression, 
ly = I1I2 ^ I4 + li- Thus, the fourth invariant is simply 

/4 = /i/2+/3 --^2detFi. (19) 

Using Ii, I2, I /a I, and I^, as obtained above with the 
knowledge of Vi , V2 , and Fi , we can apply the Simon 
separability test (Eq. (3) of our Letter). If a two- mode 
Gaussian state is entangled we know for sure that < 
and we can, therefore, fully quantify its entanglement 
either via the entanglement of formation or the negativ- 
ity/logarithmic negativity, as discussed in our Letter. 

Finally, we must emphasize that the main result of our 
Letter remains unchanged: it is possible to completely 
characterize via local operations and classical communi- 
cation (LOCC) the entanglement content of an arbitrary 
two-mode Gaussian state. Furthermore, all the results 
presented in the Letter remain valid. 

We want to thank Yang Yang, Fu-Li Li and Hong-Rong 
Li for also calling our attention on the problems related 
to Eq. (4) of our Letter while this erratum was being 
formulated. 



